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Abstract 

A calculation is presented of the plasmon and photon damping rates in a 

dense nonrelativistic plasma at zero temperature, following the resummation 

program of Braaten-Pisarski. At small soft momentum k, the damping is 

dominated by 3 —* 2 scattering processes corresponding to double longitudinal 

Landau damping. The dampings are proportional to (a/vp) 3 ^ 2 k 2 /m, where 

vf is the Fermi velocity. 
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Over the last few years, there has emerged a systematic perturbation scheme developed 
mainly by Braaten and Pisarski |I] and by Frenkel and Taylor |2| , which has finally allowed 
the calculation of some thermodynamic and dissipative properties of quasiparticles in ultra- 
relativistic plasmas (for a recent overview, see ||). It is based on a resummation method 
in which the distinction between soft scales of order gT or smaller and hard scales of order 
T is completely crucial. Basically, when the quantities to be computed depend on soft 
scales, the resummed perturbative expansion makes use of effective propagators and vertices 
accounting for long-distance medium effects as such Debye screening and Landau damping. 
The independence on the choice of gauge is guaranteed by the general properties of the hard 
thermal loops entering into the effective propagators and vertices 

One of the greatest success of these calculations has been the consistent calculation of the 
gluon damping rate at zero momentum [|J, thus solving the plasmon puzzle in perturbative 
QCD at high temperature. Also, the fermion damping rate for an excitation at rest has 
been computed || . Other interesting quantities as the damping rate of hard photons in the 
quark-gluon plasma |J, and the lifetime of a moving fermion in QED plasmas also have been 
considered within this framework 0. Only very recently, it has been shown by Blaizot and 
Iancu that the infrared divergences which plague this last quantity can be eliminated by a 
resummation based on the Bloch-Nordsieck mechanism at finite temperature || . All of these 
developments refer to a ultrarelativistic, high-temperature, zero-density regime. However, 
although degenerate plasmas (relativistic or not) are relevant in a wide area of applications 
ranging from metals to heavy ion collisions, relatively little work has been devoted to the 
computation of the damping of their excitations within the new developments [§,[L(J. In 
both of these references, the authors compute the damping of hard relativistic modes of 
momentum k ~ kp, where kp is the Fermi momentum. 

Also, it is interesting to consider the problem of the damping rate of soft modes of 
momentum and energy of order gkp and gEp respectively. They correspond to collective 
excitations in the plasma. Actually, the perturbative calculation of damping rates of these 
quasiparticles in the nonrelativistic regime is an old question. Over the sixties, there was 
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some controversy concerning the damping rate of plasmons. A preliminar calculation was 



carried out by DuBois JlT] who neglected screening effects. Ninham, Powell and Swanson 
fl"2|j repeated the calculation by including screening but missing out effective vertices which 
contribute to the same order in e as the tree-level diagrams. Lately, DuBois and Kivelson 



corrected that procedure by taking into account the screening and the effective vertices 



and this seemed to be the final perturbative answer. However, the consistency of their 
method is not guaranteed because the mixing of the soft scales (gEp, gkp) and the hard 
scales (Ep, kp) which is present in their calculations. Moreover, these authors do not report 
the numerical values of dampings. 

In this paper, I will largely follow the resummation program of Braaten and Pisarski to 
explicitly compute the plasmon and the photon damping rates at nearly zero momentum 
in a non relativistic degenerate electron gas to leading order in the electromagnetic fine 
structure constant a = gl/^K — 1/137. In this context, it is usual to introduce another 
coupling constant r s = (a/vp)(97i/4) 1 ' 3 . The degenerate limit is then r s — > 0. Resummation 



is required in this limit because the plasma frequency il p = g e kFyVF/3/ir is soft versus the 
Fermi energy Ep = k F /2m, their ratio, of order r l J 2 , playing the role of the coupling constant 
g in the Braaten-Pisarski terminology. 

The main results of this paper are the formulae for dampings rates listed below in 
Eqs. (E20|)-(f22|). The calculation, which is exact at leading order in r s , does not include 
the effects of current-density and current- current interactions which lie outside the non rel- 
ativistic approximation and are suppresed by powers of the Fermi velocity over c. However, 
it is would be noted that the non-leading terms can be computed also within the framework 
described below. 

Let us proceed to outline our calculations. The scattering processes in the plasma con- 
tributing to the plasmon and photon damping rate to lowest order in r s can be identified 
by cutting the one-loop self-energy graphs with effective propagators and vertices drawn in 
Fig. [IJ. There are three types of contributions depending on which piece of the imaginary 
part of effective propagators is picked. The pole-pole term corresponds to the decay of the 



plasmon or photon in a pair of plasmons which is forbidden by kinematics. The pole-cut 
term describes 2^2 contributions which are zero since k° = Q p is at or below threshold 
for these processes. The only nonzero contribution is the cut-cut. It comes from 3^2 
processes: j*(k) e(pi) eip-z) — > e(pa)e(p4), two electrons from the Fermi sea end up above 
the Fermi surface through absorption of the plasmon. The damping arises enterely from 
the double Landau damping in the longitudinal exchange between the electrons outside and 
inside the Fermi sea as can be seen in Fig. ^|. It should be noted that the first four scatter- 
ing diagrams arise from cuts in both effective three-vertices and effective four-vertex. The 
remaining five diagrams obtained by exchange of particles P3 <-> P4 do not enter to lowest 
order in r s , since they come from cuts of higher-loop self-energy graphs. 
The damping rate Ti t {k) corresponding to this process is 

r n\ 1 1 f ^P 1 d3 P 2 d3 Pz d3 P* 1 \ i Hi I w n l w 

VlM = 2Mfc)2!2! J (2^)W (2^ (2^^ (Pl) ^ fe)[1 " MPs)][1 " Uf ^ )] 

x (27r) 4 5(u ltt (k) + E x + E 2 - E 3 - E 4 )5^(k + Pl + p 2 - p 3 - p 4 ) , (1) 

where |-M; 5 t| 2 denotes the scattering amplitude squared, averaged over the helicity states of 
the incoming photon and summed over the spin states of the initial and final electrons. A 
double factor 1/2! is included corresponding to the two pairs of identical particles present in 
the initial and final states. The electron energies are Ei = pj/2m and the Fermi distribution 
is np(pi) = 9 (hp — Pi)- Now it is easy to write the matrix element from the diagrams in 
Fig.2, 

M l = g 3 e [g 2 e U^(P 3 - Pt, P 4 - P2)D{u, p 3 - pi)D(o/, p 4 - P2) + Jo\e (k) , (2) 
Mt(X) = gM^f\Ps ~ Pi, P± ~ P2)D{u, p 3 - Pi)D(u', p 4 - p a ) + J^k, A) . (3) 

The alert reader may wonder about the different orders of perturbation theory which seem to 
be mixed in the previous equations. In fact, we shall see below that the scattering amplitudes 
are of order g\ due to cancellations between the effective vertices and propagators. 
In the Coulomb gauge, the longitudinal effective propagator is 



where Q\ is a Legendre function of the second kind 

1 + x 



Qi(x) 



-l + -ln 
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\x\). 



(5) 



This is not the retarded propagator but coincides with it for uj > 0. Since the Fermi 
numbers in Eq. ([!]) and the energy conservation condition uj^ t {k) = uj + uj' ensure that 
energy transfer variables uj = E 3 — E\ and uj' = E 4 — E 2 are both positive, the effective 
vertices appearing in the amplitudes can be computed from their euclidean counterparts 
depending on imaginary frequencies. The complete calculation of the EE's is rather involved, 
but fortunately only the 'hard dense loop' part of them is required. In this approximation, 
one considers uj,uj',p, k -C E F , k F . Besides, the leading contribution to the ITs comes from 
a region in phase space where all internal fermionic momenta are in a narrow shell around 
the Fermi sphere, so that one can make an expansion of the integrand around the Fermi 
momentum after the imaginary frequency integration is performed. Then, the analytic 
continuation with the retarded prescription from the imaginary frequencies to uj + ie and 



uj + ie gives 
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and 



Uf\uj,uj',p,k) 



k F pj 

Tl 2 UJ t 



-Qi(uj/v F p) + Qx(uj'/v F p)) 
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uJ-^Qi(uj/v F p) +uj'-^jQ 1 (uj'/v F p) 
[ujQ^uj/vpp) + uj'Q 1 (uj'/v F )} + 0{k 2 ) 



(7) 



where p = 1/2 (p 3 — pi + p 2 — P4). Terms proportional to kj do not appear in Eq. (|7|) 



since they do not contribute because the polarization vectors satisfy the Coulomb gauge 
constraint k ■ e(A) = 0. 

The terms j7o and Jj in the amplitudes are easily written from the real part of the non- 
relativistic electron propagator because the virtual fermion is not on shell. They are of the 
form 

Jo = D(uj, p 3 - pi)— — — 3—— — — + three similar terms , (8) 

E(p 2 ) + 0Ji - E(p 2 + k) 

rr ^, \ (2pi + k j )(2m)- 1 , . 

jj = D(LJ,p 3 — pi)— r —, rr + three similar terms . (9) 

E(P2) +u t - E(p 2 + k) 

Putting this together, we find the amplitudes to lowest order in k 

~k 2 p 2 - 4k ■ p 2 



Mi = g, 



mu)f 



-D(u,p)D(u',p) 



M t (X) = g 



k ■ p (2 k • pmu' + 2p 2 ■ kp 2 - k ■ pp 2 ) ^ . , 

+ 2-3-3 D(u,p) 

m z p z uJi 

k ■ p (2k • pmu - 2pi • kp 2 - k ■ pp 2 ) , , 
+ 3-2-3 D(u ,p) 

m z p z ujf 

Ak-p 2 p { , n , , > 

D(u,p)D(u ,p) 

muj t 

2 k • p rmo' p l + k • p pi — k • p p l + k • p 2 p % 

+ 2 D(u,p) 

2 k ■ p moj p % — k ■ p p\ — k ■ p p % — k ■ pip % . 
+ 2 D ( u >P) 



e (fc), (10) 



ei(k,A). (11) 



Because of simplifications between Legendre functions, these expressions are far simpler that 
one might have expected given the complicated form for the effective vertices. It should be 
noted the cancellation of terms of first and zeroth order in k in M.\ and M. t respectively. 
Also, they agree with those previously computed, long time ago, in Ref. [f[3|| . The terms 
proportional to k ■ pp 2 and k • pp 1 can be dropped because they are subleading versus 
P1.2 ■ kp 2 and k ■ p\ )2 p l respectively. The reason is that both pi, p 2 and the final momenta, 
which for k — > 0, can be approximated for pi + p and p 2 — p, are close to the Fermi surface, 
which means that p kp. 

We still have to give the expressions for the polarization terms, 
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eo(k) = ^v^), (12) 



e j {k,\) = yJZ t (k)e j (k,\), (13) 

where Zi(k) and Z t (k) define the residue functions at the poles in the external effective 
propagators 0] and ey(k, A) are orthogonal to k and normalized so that e(k, A)-e(k, A)* = 1. 



To lowest order needed, these are Zi = Z t + 0(k 2 ). The sum over transverse polatization 
vectors is 

e;(k, A)e i (k, A)* = <J y - k i k j . (14) 

Now, it is convenient to make use of the fact that T^ t are independent of the direction 
of k. We can therefore average the squares of A4ij over the directions of k. Then, the pi 
and P2 integrals can be done using the 'hard dense loop' results 

/I* Mpi)(1 _ Mpi + p)) ,(._Ei_P_|l) , (15) 

valid for p <C kp. Lastly, the spin and polarization summations give 

x [23 wu/p 4 flj | D (w , p) | 2 1 D (a;', p) \ 2 
+±uu'{p 2 v 2 F -u 2 )\D{u',p)\ 2 

+Auu'(p 2 v 2 F -oj' 2 )\D(oj,p)\ 2 ], (18) 



T t(k) = /_ dp j o " duj^' du'5(Q p -uj-uj') 



roc 



x[16ww / p 4 fiJ|iP(w,p)| a |D(a/,p)| 2 
+3^c/(p 2 4-cu 2 )|,D(c/,p)| 2 

+3^'(p 2 4-cu /2 )| J D(cu,p)| 2 ]. (19) 

The integrals can be computed through the introduction of the dimensionless variables 
x = uj/vfP, x = uj'/vfP and q — vpq/flp. The final results are 
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Fi,t( k ) = a ltt 



( 



a \ 3 /2 k 2 



(20) 



UjJ 200m 



where the constants are determined by numerical integration to be 



~ 6.233 87, 



(21) 



at ^ 4.636 08. 



(22) 



In conclusion, I have computed the leading contribution in v F /c to the damping rates 
of soft quasiparticles in a degenerate non relativistic plasma for small momentum k. This 
contribution comes from density-density interactions in the plasma and vanishes as k — > 
0. For k — 0, there are also contributions of higher orders in Vf/c coming from velocity 
dependent interactions. These could be computed within the same framework by including 
one transverse effective interaction which would give the subleading correction. 
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FIGURES 

FIG. 1. One-loop self-energy diagrams and effective vertices. 
FIG. 2. Direct scattering graphs for 7* (A;) e(pi) e(j?2) —* e(j>3) e(P4) • 



9 




FIG. 1. One-loop self-energy diagrams and effective vertices. 




FIG. 2. Direct scattering graphs for J* (k) e(pi) e(p2) - ^ e(P3) e {Pi) 



